#! /usr/bin/env Rscript
setwd("~/Documents/facultad/tesis/goldStandard/estadisticas/")
# Estos valores se calculan con GoldStandardVsSystem del proyecto ocovinu en java.
#activacionGoldStandar = c(1.8,1.2,3.0,2.6,2.4,2.2,2.6,1.8,2.8,3.0,2.4,1.4,2.2,2.4,2.0,2.6,2.6,2.0,1.0,1.6,2.4,1.2,2.2,2.2,2.8,2.0,2.4,2.2,2.8,2.4,2.6,2.2,1.8,1.2,2.8,2.0,2.4,2.6,2.0,2.4)
#agradabilidadSistema = c(2.1,2.05,2.36,2.55,2.37,2.3,2.29,2.37,2.26,2.5,2.56,2.42,2.22,2.32,2.49,2.42,2.49,2.43,2.18,2.52,2.15,2.3,2.14,2.2,2.35,2.25,2.12,2.38,2.35,2.31,2.17,2.32,2.33,2.43,2.22,2.29,2.48,2.33,2.26,2.26)
#imaginabilidadGoldStandar = c(1.8,1.2,3.0,2.6,2.4,2.2,2.6,1.8,2.8,3.0,2.4,1.4,2.2,2.4,2.0,2.6,2.6,2.0,1.0,1.6,2.4,1.2,2.2,2.2,2.8,2.0,2.4,2.2,2.8,2.4,2.6,2.2,1.8,1.2,2.8,2.0,2.4,2.6,2.0,2.4)
#agradabilidadGoldStandar = c(1.8,2.0,2.6,2.8,2.8,1.8,2.4,1.8,2.6,2.4,3.0,2.0,1.2,1.8,2.8,2.2,2.2,2.4,2.0,2.2,1.4,1.6,1.0,1.2,2.2,2.4,1.4,2.4,2.4,1.8,2.2,2.2,1.8,2.4,2.8,2.4,2.8,2.4,1.2,2.2)
#activacionSistema = c(2.1,2.25,2.39,2.19,2.38,2.5,2.08,2.31,2.35,2.5,2.36,2.54,2.34,2.26,2.29,2.25,2.42,1.93,2.0,2.12,2.37,2.13,2.08,2.55,2.37,2.19,2.21,2.34,2.47,2.23,2.19,2.23,2.27,2.42,2.35,2.19,2.47,2.38,2.22,2.43)
#imaginabilidadSistema = c(2.83,2.25,2.72,2.62,2.71,2.53,2.69,2.45,2.26,2.68,2.64,2.51,2.44,2.54,2.7,2.6,2.63,2.66,2.7,2.53,2.5,2.71,2.53,2.67,2.45,2.59,2.26,2.27,2.69,2.72,2.51,2.54,2.31,2.42,2.41,2.56,2.83,2.6,2.65,2.54)
textos = c(1:40)
agradabilidadGoldStandar = c(3.0,2.8,2.8,2.8,2.8,2.8,2.6,2.6,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.2,2.2,2.2,2.2,2.2,2.2,2.2,2.0,2.0,2.0,1.8,1.8,1.8,1.8,1.8,1.8,1.6,1.4,1.4,1.2,1.2,1.2,1.0)
agradabilidadSistema = c(2.56,2.55,2.37,2.49,2.22,2.48,2.36,2.26,2.29,2.5,2.43,2.25,2.38,2.35,2.43,2.29,2.33,2.42,2.49,2.52,2.35,2.17,2.32,2.26,2.05,2.42,2.18,2.1,2.3,2.37,2.32,2.31,2.33,2.3,2.15,2.12,2.22,2.2,2.26,2.14)
activacionGoldStandar = c(3.0,3.0,2.8,2.8,2.8,2.8,2.6,2.6,2.6,2.6,2.6,2.6,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.2,2.2,2.2,2.2,2.2,2.2,2.0,2.0,2.0,2.0,2.0,1.8,1.8,1.8,1.6,1.4,1.2,1.2,1.2,1.0)
activacionSistema = c(2.39,2.5,2.35,2.37,2.47,2.35,2.19,2.08,2.25,2.42,2.19,2.38,2.38,2.36,2.26,2.37,2.21,2.23,2.47,2.43,2.5,2.34,2.08,2.55,2.34,2.23,2.29,1.93,2.19,2.19,2.22,2.1,2.31,2.27,2.12,2.54,2.25,2.13,2.42,2.0)
imaginabilidadGoldStandar = c(3.0,3.0,2.8,2.8,2.8,2.8,2.6,2.6,2.6,2.6,2.6,2.6,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.4,2.2,2.2,2.2,2.2,2.2,2.2,2.0,2.0,2.0,2.0,2.0,1.8,1.8,1.8,1.6,1.4,1.2,1.2,1.2,1.0)
imaginabilidadSistema = c(2.72,2.68,2.26,2.45,2.69,2.41,2.62,2.69,2.6,2.63,2.51,2.6,2.71,2.64,2.54,2.5,2.26,2.72,2.83,2.54,2.53,2.44,2.53,2.67,2.27,2.54,2.7,2.66,2.59,2.56,2.65,2.83,2.45,2.31,2.53,2.51,2.25,2.71,2.42,2.7)

png(file="sistemaVsGoldStandard-Agrado.png", width=450, height=400)
plot(agradabilidadSistema, type="o", col="blue", ylim=c(1,3), xlab="Textos", ylab="Agrado")
lines(agradabilidadGoldStandar, type="o", pch=22, col="red")
abline(lm(agradabilidadSistema~textos), col="blue")
abline(lm(agradabilidadGoldStandar~textos), col="red")
title(main="Sistema vs Gold Standard (Agrado)", col.main="black", font.main=4)
legend("topright", c("Sistema","Gold Standard"), cex=0.8, col=c("blue","red"), pch=21:22)
dev.off()

png(file="sistemaVsGoldStandard-Activacion.png", width=450, height=400)
plot(activacionSistema, type="o", col="blue", ylim=c(1,3), xlab="Textos", ylab="Activación")
lines(activacionGoldStandar, type="o", pch=22, col="red")
abline(lm(activacionSistema~textos), col="blue")
abline(lm(activacionGoldStandar~textos), col="red")
title(main="Sistema vs Gold Standard (Activación)", col.main="black", font.main=4)
legend("topright", c("Sistema","Gold Standard"), cex=0.8, col=c("blue","red"), pch=21:22)
dev.off()

png(file="sistemaVsGoldStandard-Imaginabilidad.png", width=450, height=400)
plot(imaginabilidadSistema, type="o", col="blue", ylim=c(1,3), xlab="Textos", ylab="Imaginabilidad")
lines(imaginabilidadGoldStandar, type="o", pch=22, col="red")
abline(lm(imaginabilidadSistema~textos), col="blue")
abline(lm(imaginabilidadGoldStandar~textos),col="red")
title(main="Sistema vs Gold Standard (Imaginabilidad)", col.main="black", font.main=4)
legend("topright", c("Sistema","Gold Standard"), cex=0.8, col=c("blue","red"), pch=21:22)
dev.off()

cor.test(agradabilidadGoldStandar, agradabilidadSistema)
cor.test(activacionGoldStandar, agradabilidadSistema)
cor.test(imaginabilidadGoldStandar, agradabilidadSistema)
cor.test(activacionGoldStandar, activacionSistema)
cor.test(imaginabilidadGoldStandar, activacionSistema)
cor.test(imaginabilidadGoldStandar, imaginabilidadSistema)

